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New hypergeometric formulae to tt arising from M. Roberts 

hyperelliptic reductions 

Giovanni Mingari Scarpello * Daniele Ritelli ^ 


Abstract 

In this article we developed a special topic of our pure-mathematics papers mHEiii] concerning the 
hypergeometric theory. Based upon a Roberts’s reduction approach of hyperelliptic integrals to elliptic ones, 
[5], and on the simultaneous multivariable hypergeometric series evaluation of them, several identities have 
been obtained expressing tt in terms of special values of elliptic, hypergeometric and Gamma functions. By 
them TT can be provided through either only one or two parameters and through the imaginary unit. 

In any case, such results, all unpublished and undoubtably new, will provide, beyond their own beauty, 
a useful tool in order to check the routines (more or less naive) which one can build for the practical 
computations of Lauricella’s functions met frequently in researches on Mechanics or Elasticity, for instance 
by ourselves in [SlIzlIllE]. 


1 Introduction 

First of all, a short framework of the Roberts’s booklet, [S], which is the start-up of all our next processing. 
After, a previous hypergeometric set of statements will be provided. 


1.1 A bit of history 

The british mathematician Michael Roberts (1817-1882) was attired by the hyperelliptic integrals since the 
publication of some relevant post-Jacobi papers by Riemann and Weierstrass. In 1871 he published a short 
tract, [S], which collected his lectures on the subject. He constructed a trigonometry of hyperelliptic functions 
analogous to the well-known one for the elliptic functions. Its sixth chapter deals with the first class hyperelliptic 
integrals depending on, or which can be reduced to, elliptic ones. They include for instance some integrals like: 


I 


\/Qg{x) 


dx 


( 1 . 1 ) 


where Qg{x) is a square-free polynomial of any given degree 5 > 4. The allowable power k has to be determined 
by analysis of the possible pole at the point at infinity on the corresponding hyperelliptic curve. When this is 
done, one finds the condition to he k < g — 1, or, in other words, k at most 1 for degree 5 or 6 of Qg{x); at 
most 2 for degrees 7 or 8, and so on. 

Next Roberts gives the analogous with Fagnano’s theorems founding upon the analogy between conic sections 
and the lines of curvature of an ellipsoid, linking to his best loved research subject of geodesics on such a surface. 
A commemorative article, m, issued two years after his death, describes the scientific profile of M. Roberts. 

In his booklet [5], section 63, Roberts takes into account reductions to elliptic of some special hyperelliptic 
integrals as: 


Rn = 


\/x® — px^ + qx^ — px^ + 1 


dx 


( 1 . 2 ) 


where n takes the integral values 0, 2, 4. Reduction is carried out through the variable transformation: 


/(x) := X -I— = u 
X 


X = 


- (^u± _ 4 ^ 


dx = 


1 ± 


— 4 


(1.3) 
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The double sign is depending on being f{x) not monotonic and + has to be taken when x G [0,1] and — if 
X G (1, oo). 


U 



Figure 1: Cauchy-Schlomilch variable transformation 


If the integration domain covers both the intervals [0,1] and [l,oo), integration has to be properly split. 
After (1.3) one finds: 


Rn = 


1 f (\/u2 — 4 ± u) 


n —1 


1 


— {p + 4:)u‘^ + q + 2p + 2 


du. 


(1.4) 


Such a transformation used by Roberts is known as Cauchy-Schlomilch one, useful to evaluate many non¬ 
elementary definite integrals, most of which cannot be obtained by symbolic packages or otherwise. In El the 
(1.31 is a tool for computing several of them and some historical notes are provided. In [T] (1.3) is also extensively 
employed, following its use introduced by A.M. Legendre. In order to perform the necessary reduction, after the 
powering, one gets two integrals. The first holds the radical _ 4 - j^^s also to be transformed by putting 
V = Restricting for shortness to the case x G [0,1], the Roberts reductions are0 



dtt 


du 


y/u'^ - (p + 4:)u‘^ + q + 2p + 2 

du 

y/u^ - {p + 4:)u'^ + q + 2p +2 

(1 — u^)du 


— {p — 4:)v^ + q — 2p + 2 j 
I" du 

J y/v'^ — {p — 4:)v'^ + q — 2p + 2 ^ 
(l-l-u^)du \ 


y/u"^ — {p + 4)u2 _(_ g _|_ 2p -I- 2 J — {p — 4)u^ -I- g — 2p -I- 2 


(1.5) 

( 1 . 6 ) 
(1.7) 


where the left hand side hyperelliptic integrals have been reduced to a combination of elliptic ones with quartics 
under the quadratic radical. We assumed to move to the expressions (1.5), ([13 because in [ 0 , 1 ] the 

transformation (1.3) is decreasing, so that, when passing to definite integrals, we will change the integration 


limits. Of course, when integrating in the interval [1, 00 ) the sign will be properly fit. Through our method of 
double evaluation we already employed in El 13 Hill], we will establish further links among tt, some multivariable 
hypergeometric functions and elliptic integrals. 


In order to treat hypergeometrically the hyperelliptic integral at left hand side of (1.2), it will be necessary 
to solve the eighth-degree equation: 

P{x) = x® — px^ + qx^ — px^ -|- 1 = 0. (1.8) 


By the aforementioned (1.3) we get: 


+ 9 ) =a;4((u"-2)"-2-p(u2_2)+< 7 ) ((^ 3 ) 


^Notice that he is not concerned with defining exactly how the integration limits are transformed and neglects integration limits 
with the integrand > 1. 
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where, due to (1.3), we put: 


1 


X -\— = u 

X 


4 = »" - 2 


+ ^,={u^-2)^-2 


Then the roots of (1.8) are found by solving the auxiliary equation: 

u'^ - {p + 4)u^ + 2p + q + 2 = 0, 


(1.8:) 


and then one goes back to originary variable x via the (1.3). Due to the symmetry of equation (1.8), which is a 


reciprocal one, if it has the root r, then —r e 1/r will necessary be roots too. Three main cases can then occur: 

• 8 real roots • 4 pairs of conjugate complex roots 

• 4 real and 2 pairs of conjugate complex roots 


Furthermore, if 1 is a root of (1.8) it shall be double. 


In our article, by comparing the hypergeometric approach for evaluating i?„, see lemma [2.1[ lemma [272| 
and lemma |2.3[ versus the Roberts method, we will get not only several reductions of hyperelliptic integrals 
to elliptic ones, but also new formulae for tt to be added to those we presented in [1 [3 E H. The treatment 
will depend on the index n of the integral for n = 0, 2, 4 the relevant complete integrals are studied 

taking as integration limits the roots of ( |1.8[ ), or, in the complex case, by integrating on the straight line [0, oo). 
Roberts’s reductions are possible also for integrals with even further exponents and which can too be reduced 
to elliptic ones. But they have been omitted both for shortness and for avoid analogous developments but with 
much involved formulae. For instance, taking n = 6 analogously to (|1.5[), (1.6) and (|1.7|) one arrives at: 


Re — 


(—1 + 3u^ — u^)dr 


(1 + 3v‘^ + ?;^)dt 


— (p + 4)u^ + q + 2p + 2 J \/v'^ — [p — ^)v‘^ + q — 2p-\-2 


(1.9) 


where the integrals at right hand side can be given as a combination of elliptic ones. 


2 Hypergeometric forewords 

In such a section some integration formulae will be established which make use of multivariable hypergeometric 
functions and to be compared with Roberts reductions given in [ 3 . 

By treating hypergeometrically the integrals like: 


a/ — a?){x'^ — b'^){x‘^ — — (P) 


dx 


( 2 . 1 ) 


several assumptions can be made: the numbers a, 5, c, d can be real and 0<a<b<c<doTa, b real with 
0 < a < b and c = d or a, b, c, d complex with a = b, c = d. When necessary, suitable conditons on s will 


be required so that the integral (2.1) succeeds to be convergent. Our first lemma 2.1 deals with a, 6, c, d all 
real, and further ones and |2.3| will deal with complex roots. The special functions necessary in order to 
evaluate the integrals (2.1) are the hypergeometric Lauricella functions of n G N’*' variables, see 
m. defined as: 


and 


D 


a, 6 i,...,&„ 


Xl, 


■= E 


(a) 


miH- \-'mn ipl)mi 


(^n)r 


mi . ,m„ GN 


(c)mi+-+m„TOi! • • -mnl 


where {x)k is the Pochhammer symbol, and under the hypergeometric series usual convergence requirements 
|xi| < 1,..., \xn\ <1. If Rec > Rea > 0, the relevant integral representation theorem provides: 


D 


a, 61 ,..., 


X\ , . . . , Xyi 1 — 


r(c) 


,a—1 


(l-u) 


c—a—1 


dw 


r(a) F(c — a) Jq (1 — xiu)^^ • • • (1 — 


( 2 . 2 ) 


allowing the analytic continuation to C" deprived of the cartesian n-dimensional product of the interval (1, 00 ) 
with itself. 

In order to define the hypergeometric integration, we take apart the different cases of roots of the eighth 


degree equation under the radical in ( 2 . 1 ). 
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2.1 The case of 8 real roots 


will see the following basic definite integrals: 

j[''\a,b,c,d) = [ - 
JO \ 

J^''\a,b,c,d) = 
j'^^\a,b,c,d) = 


'b \J (a ;2 — a^)(a ;2 — b'^){x^ — c'^){x‘^ — d?) 

_ 

'd yj (a:^ — a^){a :2 — b‘^){x‘^ — c^){x‘^ — d?) 


- 6 ^)(x^ — c^)(x^ 

- d^) > 0. We 

—™ dx 

(2.3) 

d?) 

-dx 

(2.4) 

d?) 

-dx 

(2.5) 

- d^) 


By the next lemma 2.1 we will be allowed to compute ( |2.3| ), (2.4) and (2.5) hypergeometrically: 
Lemma 2.1. If0<a<6<c<ci, then: 


/^Habcd)- 

(a,b,c,dj- 1 ^ 

J^"‘\a,b,c,d) = 

J^"\a,b,c,d) = 


1+s . Ill 
2 ’ 2 ’ 2 ’ 2 


2 2 2 
a a 


2^(62 - a2) (d2 - 62) 


I s- 1-2 

\ 2 

/ 1 ■ iizs 1 1 

^(3) / 2 ’ 2 ’ 2 ’ 2 

1 


62 ’ c2 ’ d2 I 

,2 U2 


_ a/^L (^) ^(3) / 2!2’2>2 

2 ^(d 2 _ a 2 ) (d 2 _ 52 ) (^2 _ c 2 ) r( 2 -f) ^ I 2 


62 ’ 62 — a2 ’ (62 _ 52 I 

2 62 


(2.3d) 


Mb) 


a 


— cP^ — d? 




where in (2.3) and (2.3D)we assume s > —1 and in (2.5) and (2.5 dI we assume s < 3. 

Proof. Through the integral representation theorem and our customary approach as in papers Emuiii] , we 
will compute the integrals (2.3), (2.4), (2.5). Formula (2.3d) can be obtained via the change x = Uy/u in 
(|2.3|) and calling (2.2). Analogously (2.4dI is found putting x = yjb'^ + (c2 — b'^)u, whilst (2.5|d|) follows from 

X = d/yj\ — U. □ 

2.2 The case of 4 real roots and 2 pairs of conjugate complex roots 

The following lemma will concern the integrand with complex roots. 

Lemma 2.2. Let it be a, 6 G R, a < 6, z G C \ R. We will consider the integrals: 


Mi(a, 6 , z) = 
M 2 (a, 6 , z) = 


'0 yj (a ;2 — o 2 )(a ;2 — b‘^){x'^ — z‘^){x'^ — z^) 
r I" 


'b yj{x^ - a2)(x2 - 62)(cc2 - z2)(a;2 - z^) 


dx 

dx 


where in (|2.6|) we assume s > —1 amd in (|2.7|) s < 3. Then: 

of of 


Mi(a, 6 , z) = 
M 2 (a, 6 , z) = 


r (^) (3) 

r T 


s+ 1 . Ill 
2 ’ 2 ’ 2 ’ 2 


s-|-2 

~T~ 


26|zM(^) ^ 

__ r (I - f) _^( 3 ) j 2 > 2 ’ 2 ’ 2 

2 v /(62 - a 2)(62 - z 2)(62 - z ^) T (2 - f ) ^ 


1.111 


4—s 

2 


( 2 . 6 ) 

(2.7) 


Z 2 


@) 


Proof. On the steps of proof of lemma 2.1 in the integral (3.22i) the change x = Oy/u is used, whilst for (2.7) 
the change x = 6(1 — u)“^/^. Thesis (2.6d) and (2.7 d) follow from the integral representation theorem. □ 
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2.3 The case of 4 pairs of conjugate complex roots 

Finally, analogous of two preceding lemmas [2.1[ |2.2| when the integrand has all its roots complex: 
Lemma 2.3. Let a, 6 G C \ K then, if s < 3 


L{a, b) = 


° \J — b^) 


dx 


^ i i i i 

7r(,i — S) (4) / 2 ’ 2’ 2’ 2’ 2 


4cos(f) 


7ZS ^ D 


( 2 . 8 ) 


1 + a^, 1 + a^, 1 + 6^, 1 + 6 


Proof. In order to apply (2.2), it will be enough a change in (2.8) putting first x = -^pv and after v = (\ — ufju 


in such a way obtaining the integral on [0,1]: 


u 2 (1 — u)- 


2 7o a/I — (1 + a^) upl — (1 + bp upl — (1 + c^) upl — (1 + d '^): 


: du. 


□ 


3 The case n = 0 


We want to study the consequences of assuming the values n = 0, 2,4 in (1.21 and for each of them the relevant 


sub-cases. At each n-value will be devoted a separate section. We start with n = 0. 


3.1 Sub-case: 8 real roots 

In order to use the Roberts transformation, due to the P{x) structure, we will have two numbers a, b such that: 


P{x) = - a^)(x^ - ^){x‘^ - bPix^ - 




Then in the (1.2) integrand, we have: 


ia^ + bp il + Pbp 1 + PP + (a2 + bp ^ 

P = —^2 -^ 


(3.1) 


(3.2) 


Fixing two real numbers 1 < a < 6, by we get complete elliptic integrals to be integrated in [0,1/&], [1/a, a], 
[b, oo), when P{x) has the form Let us see them apart. 

Interval [0,1/b] 


We introduce a special symbol for describing how (1.2) specializes taking into account: the index n, the inte¬ 


gration interval, and the number of P(x)’s real roots. We put: 




x"dx 


'o ^{x^-aP{x^-P){x^-bP{x^-P) 


(3.3) 


The left index “8” recalls the number of real roots of (1.8); the apex concerns the integration interval; the right 


the generators of the roots of (1.8). If n = 0 we have: 


index marks the degree of power at numerator of (1.2), and finally the argument within round brackets defines 


Theorem 3.1. With 1 < a < 5 we get: 






fb{l + ap 




fb{l-ap 


1 

2 \l + b^“ \a{l + bPJ l-b^“\a{l-bP 


(3.4) 
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Proof. Let us compute the last above integral by the Roberts reduction formula (1.5), through (3.2). We get: 




dw 


V 


'b+i 


+ 


du 




{b+lf -u^ 


[a — - ) — 




Both integrals at right side of (3.5) are complete elliptic of first kind: 

r dx 


/ 

(3.5) 

(3.6) 

□ 

By means of lemma |2.1| we can get our further formula providing tt and to be added to other ours given in 

I1I3II1I1]. 




a 


Ja a/ (a^ — x2)(/32 — x'^) Oi 
see entries 3.151-12 p. 277 of [2] and 215.00 of [IS]. Thus thesis follows. 


, P <01, 


Theorem 3.2. If 1 < a < 6 then 


K 


(a^-l-l)6\ 5^ 


-b 


K 


{a^-l)b" 


b"^ + 1 \ {b"^ + l)a I b'^ — 1 I (6^ — 1) 


f(3) / 

1.111 

2 ’ 2 ’ 2 ’ 2 


1 1) 

" 1 

1 

62’ 

q2^2 ’ Ij 4 j 


(3.7) 


Proof. We use the first integration formula of lemma 2.1 for a integral like (2.3). The roots identihcation is: 


a 

6 

c 

d 

t 

1 

1 

t 

lib 

1/a 

a 

b 

1 true 

for 1 

< 

a < b 


bR 


[ 0 , 1 / 6 ] 

0 


1.111 




1 1 


’ q, 2^2 ’ fji 


(3.8) 


Thesis (3.7) follows by reducing the right hand side of (3.5), applying the (3.8) for computing 8i?o°’^^^^(a, 6) 
hypergeometrically and then getting tt. □ 

When one of P{x) real roots has the value 1, such a root has to be double according to the order reduction 

(3) 

between the three variable hypergeometric function and the two variable one named Fi: 


/ a; 61 , 62 , 63 


D 


x,x,z\ = Fi 


a; 61 -I- 62, 63 


X, z 


and also due to the relationship: 

So that one is allowed to establish: 
Corollary 3.3. If 6 > 1 then 


K{0) = 


dx 


2 


262 


(1 + 62) 2 (62 - 1) Fi 


1-1 1 

2 ’ ’ 2 
1 




(3.9) 


(3.10) 


(3.11) 
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Interval [1 / a, a] 

Now let us consider the complete integral on [l/a,a] : 


h /( 


dx 


x2-a2)(a;2 - 1 )(x2-62)(^2_ 1 ) 


(3.12) 


(3.13) 


By the trasformation (1.3) we can establish the following integration formula: 

Theorem 3.4. Let it be 1 < a < 6 , then 

sRo K^)- ■ 

Proof. Integral (a, ^) will now be computed via transformation (1.3): for the purpose, the integration 

interval has to be split: 

d[ . 

8^0 




da; 


da; 


= /W(a,6) + /f)(a,6) 

For I^\a, b) one goes on with transformation (1.3), namely, moving there x within [0,1], the root 

a; = - — \/ 1(2 — 4^ . 

has to be chosen. Therefore, after changing the integration limits, we get: 

- I ... 


l'i\a,b) = - 


dat 


du 


V 






and through formulae 3.152-8 and 3.152-7 p. 276 of [14] or entries 220.00 and 219.00 from[15j: 


^ 2 &2 1 - (^2 + 1) (52 _ 1) > „ (^2 + 1) 


b ( ~ 1 ) (b^ + 1 ) b (of -f 1 ) 

F arcsin ■ 


6 j^(b{a^-l) 


52-1 


.(62-1) 


/n't 

For I 2 {a, b) the approach is similar: through the transformation a; -I- 1/a; = at with x within [1, a] with a > 1, 
so that the root to be chosen is: 

1 


As a consequence: 

( 


l2^\a,b) = ^ 




^=2 


dat 


(^u + — 4^ 


dat 


V 




(a+i) -aa 2 j [( 6 +i) - aa 2 j 

= 1 ( _Lcsin.is!jJiif±il +^k( 

2 62 + 1 y (a 2 + 1 ) (62 - 1 ) ’ a (62 + 1 ) y ^ 62 - 1 ^ a (62 - 1 ) ^ ^ 

■ (‘2^ (2\ 

Thesis (3.13) follows by adding the expressions obtained of /) {a,b) and I 2 '{a,b). 




□ 


As a consequence of theorem 3.4 just proved, by lemma 2.1 again, equation (2.4 d), we obtain the second tt 
formula of this paper: 
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Theorem 3.5. If 1 < a < 6 : 


2 -^— a^) { 0 ?^^ — 1 ) 


K 


a ( 6 ^ — 1 ) 


t 3 ('^2 _ 2^^ 


/ i- 1 i 1 

p(3) / 2> 2’ 2’ 2 

1 


l-a^ 


4(1 — a^)b^ — 1 


(3.14) 


b^ — ’ a^b'^ — 1 


Proof. Through lemma |2.1| equation (|2l4p| we compute the integral (|3.12| so that: 

(1 - 0^)62 - 1 




^bir 


2 \/ (&2 — a'^)(a‘^b'^ — 1) 


/ i- i i i 
t( 3 ) I 2 ’ 2 ’ 2 ’ 2 


D 


1 


1 — a"', 


b'^ — ’ a^b"^ — 1 


By comparing (3.12 d) with ( |3.13 1, thesis (3.14) follows. 
Interval [b, 00 ) 

Now let us deal with the complete integral along [ 6 , 00 ) : 


( P^ ) 

□ 


/ OO 

V 


da; 


(^2_a2)(^2_ 1)(3,2_52)(3,2_ 1) 


Theorem 3.6. If 1 < a < & then 




K 


6 (o 2 - 1 ) 


K 


2\b'^-\ \a{f3^-\)j 52 + 1 \a{b^ + l) 


&(a 2 + 1 ) 


(3.15) 


(3.16) 


Proof. By the Roberts reduction again, being the region where the transformation (1.3) is growing, we find: 

I ... 




du 


du 


\ 


[{- -ir- [{b -ir- d[{-+ir- [(^ +ir--^\j 


(3.17) 

Then thesis (3.16) follows by computing the complete elliptic integrals in (3.17) through the entries 3.153-11/12 
p. 277 of [15^7215.00 and 216.00 of [E]. □ 


From lemma 2.1 through (2.5 3 ) for computing the integral (3.15) and replacing the outcome with the 


theorem 3.6 formula for tt follows: 
Theorem 3.7. Let a < 6 < 1, then: 


2 /(62 — a2)(l — 0252) 

1 -p 




5(a2 -I- 1) 
0(52 -I- 1) 


/ 1.111 

p(3) / 2’ 2’ 2’ 2 

2 


1 1 a2 

1 — 6'^ ’ 1 — 0252 ’ Qr2 _ ^2 


(3.18) 


Proof. By lemma 2.1 equation (2.5 3 ) we have: 

dx 


poo 

, 4 ’“)(a, 6 )=/ - 

A ^ 


nab 


(^ 2 _^ 2 )(^ 2 _^ 2 )(^ 2 _^)(^ 2 _^) 

/ 1.1 1 1 


1 


U3) 2’ 2’ 2’ 2 


1 1 a2 

1 — 6 ^ ’ 1 — 0252 ’ a 2 — 52 


4 v'(64 _ 1 ) (62 - a2) (0252 _ 1 ) ^ 

Our further tt formula follows after equating to ( 3.19| ) and comparing to (3.16). 

Taking the limit for a —3 1, a formula to tt comes out characterized by only one parameter. 


(3.19) 


□ 
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Corollary 3.8. If & > 1 then 


2(62-1)' 


K 


2b 


1 + 62 


^54 _ 1 _ Fi 


1 1 
1-62’ l_fj4 


(3.20) 


3.2 Sub-case of 4 real and 2 pairs of conjugate complex roots 


Let the polynomial (1.81 have 4 real roots and 2 pairs of conjugate complex ones, assumed to belong to the 
unity disk. Then we will have with a, 6 € K, 6 > 1: 


p(^) = (^^2 _ (^2 _ ^2) (^2 _ g-2*a) (^2 _ g2*a) 

= ~ ~ (1 ~ 2 cos 2a + a;^) 


(3.21) 


= 1 - 


1 + 26^ cos 2a + b‘^ ^ r, 1 + 26^ cos 20 + 6"^ 4 1 + 26^ cos 2a + 6"^ g g 


-a;2 + 2- 


-a; — 


-x“ + a;° 


62 62 62 
Restricting to the positive half straight line, if a; S [0,1/6] U [6,00) we have P{x') > 0. We pass to compute the 
integrals: 

1 

da: 




/ OO 

;- 

- 


0 (a:2 — T) (a:2 — 62) (x2 — e“2*a) (2;2 — e2“) 

da; 


— p) (a;2 — e“2*a) (^2 — e2*“) 


(3.22) 

(3.23) 


by comparing the hypergeometric evaluation of ( 3.22[ ) and ( |3.23 ), with those obtained as complete elliptic 
integrals through the trasformation (1.31. 


Through lemma 2.2 we get two further formulae for tt, we state the following theorems 3.9 and 3.10 


Interval [0,1/6] 

Theorem 3.9. If a, 6 € K, 6 > 1 then: 


/ / 26 cos a 

T+^ 


K 


26 sin a 


Vl — 262 g.Qg 2a -(- 54 


F(3) / 

1.111 

2 ’ 2 ’ 2 ’ 2 

1 6^“^ 


1 + 62 


Vl — 262 cos 2a + 64 

" l 

1 

64’ 62 ’ 62 j 

V 





(3.24) 


Proof. Through the formula (2.6 a) we get the hypergeometric evaluation of 4i?Q°’^^^^(a, 6) put into (3.22): 


/ 1.111 

p[0.1/b]/ .N - /Lp(3) / 2>2>2>2 

4rto ~ 26 ^ ^ \ 1 


64’ 62 ’ 62 


(3.22i) 


Furthermore, by transformation (|1.3|) to (|3.22|), through arguments like those seen for theorem 3.2 we get: 


pOO 

da , r- 

\ 

du 


^((6 + ^)^ — a2^ (4 cos 2 a — a2) \ 

1 ^a2 - (6 - i)^^ (4sin^ a + u2) ^ 


(jga) 

Both elliptic integrals at right hand side of (3.223) can be evaluated by entries 3.152-12 and 3.152-6 p. 277 of 
[13] so that: 


4i?r/''(a,6) = 




26 cos a 


K 


2\l + 62 \l + 62y Vl — 262 ^Qg 2a + b^ V \/l — 262 cos 2a + 64 


26 sin a 


( 3^ ) 


Thesis (3.24) comes by equating (3.22i) with (3.22:) and then expressing tt via a, 6 and the imaginary unit. □ 
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Interval [6, oo) 

Theorem 3.10. If a, 6 S K, 5 > 1 then 

2y/b'^ — l\/l — 2 cos 2a 6^ + b'^ 


F(3) / 

1.111 

2 ’ 2 ’ 2 ’ 2 

1 1 

1 ) 

" 1 

2 

1 _ 54 ’ 1 _ ^2g-2/a ’ 1 

_ ^2g2/a 1 


' f 25 sin a 

V Vl — 2 cos 2a 5^ + 5^ 
\/l — 2 cos 2a 52 + 5^ 




25 cos a \ \ 

1 + 52 


V 


1 + 52 


/ 

(3.25) 


Proof. By the hypergeometrical side we compute the integral (|3.23| by means of lemma 2.2 eq. (|2.7|), obtaining 

biT 


4R^^’°°\a,b) = 


d3) 


1.111 
2 ’ 2 ’ 2 ’ 2 


4\/5^ — l-\/l — 2 cos 2a 52 + 5^ ^ V ^ 

On the other side, the Cauchy-Schlomilch transformation provides: 

^ / 

v[ 6 ,oo 


1 


1 


1 


1 


4i?r'-^(a,5) = 2 


du, 


1 _ 54 ’ 1 _ ^2g 2 ia ’ 2^ _ ^2g2 


du 


(3.26) 


^ ^ ^u2 — (b— ( 4 sin^ a -I- 11,2) ^(^(b + ^)^ — (4cos2 a — u^) ^ 


(3.27) 


The above elliptic integrals in (3.26[) are listed in [2] entry 3.152 - 12 p. 277 so that 


4 !?, 


[6,00) 


(a>^) = o 


25 sin a ^ 5 ^^25 cos a 

2 V t/l — 252 cos 2a -I- 5"* V Vl — 252 gQg 2a -|- 5^ / 1-1-52 (^2-1-52 


K 


Thesis (3.25) follows by equating (3.26) and (3.28) and getting tt. 


(3.28) 


□ 


Remark 3.11. Notice that it is not possible, in statements of theorems 3.9 and 3.10 to make 5 — )■ 1+ because 
both integrals (3.22) and (3.23) are divergent if 5 — >■ 1+. 

3.3 Sub-case of 4 pairs of conjugate complex roots 


Even if the polynomial (1.8) does not have any real root, we can refer to a representation formula to tt, theorem 


3.12 Let all the P(a;)’s roots stay within the unity disk. This assumption, joined to the special P{x) structure 


will imply: 

Fix) = [x^ - e2“)(x2 - e-2“)(x2 - e^^'’){x^ - e'^*^) 

= (1 — 2 cos 2a x"^ -I- x'^) (l ~ 2 cos 25 a;^ -I- a;^) (3.29) 

= 1 — 2 (cos 2a + cos 25) a:^ -I- (2 -|- 4 cos 2a cos 25) x"^ — 2 (cos 2a + cos 25) x^ + a:® 

where a, 5 are two real numbers upon which later we will do some further assumptions. 

First, again, an hyperelliptic integral computed hypergeometrically: by means of lemma |2.3| this will be done 
for the integral: 

/•oo 

(3.30) 


0^0 


[ 0 ,c 


pOO 

= / 7(T 


dx 


— 2cos2aa;2 -|- x'^) (1 — 2cos25a:2 -|- x^) 


Interval [0, cxj) 


Following Roberts, applying again the theorem 3.5 then (3.30) is promptly seen as a first kind complete elliptic 
integral: then by the usual comparison of different reckoning, we get: 

Theorem 3.12. Let a, 5 S M be such that sin a > sin 5. Then: 

^ j \J sin^ a — sin^ 5 \ 

1 V 


sm a 


2 sin a 


<^(4) 

D 


3.1111 
2 ’ 2 ’ 2 ’ 2 ’ 2 


(3.31) 


1 -b 6^“, 1 -b e-2“, 1 -b 1 -b 
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Proof. The Roberts approach, built on the Cauchy-Schlbmilch transformation (1.31, has to be applied splitting 
the integration domain [0, oo) in [0,1] and [1, oo): going on with the theorem 3.5 we find that the integral (3.30) 
is given by: 


POO 

oR^°’°°\a,b) = -r 


4 sin^ a -I- it^) (4 sin^ 6 -I- it^) 


but (3.30 dI comes by entry 3.152-1 p. 276 of [Tl], so that 


o4°’~^(a,6) = 


1 


2 sin a 


K 


\/ sin^ a — sin^ h 


sm a 


( Slot :) 


On the other side the same integral is computed also through the lemma 2.3 as L(e“,e*“) in equation (2.8 1 . 
Thesis (3.31) follows as before by comparison. □ 

Remark 3.13. If a —>■ 5 then theorem |3.12| does not provide any more a tt expression because the elliptic 
integral degenerates in an elementary one 


lim 

a —>-6 


du 


du 


0 y(4sin2a + u2) {4smH + u^) dsin^a + w^ 4sina 


whilst the Lauricella function has two couples of equal arguments so that it collapses in a two-variable 
hypergeometric functions of Appell according to the general relationship: 


q4) ( a;6i,62,&3,&4 


D 


x,x,y,y] = Fi 


a; bi b2, 64 


x,y 


But, in the specific case a = 3/2, 64 = 62 = ^3 = ^4 = 1/2, c = 2 and therefore, being the sum of both b 
parameters equating the parameter c, we thus have the further reduction: 


3/2; 1,1 


x,y \ = 


2 (yi - y - - x) 

(x - y)y/{l-x){l-y) 


which, written for a; = 1 -I- e^“, y = 1 -I- e and compared with the right hand side leads simply to the Euler 
identity: 

2 1 


\/—e ^ -|- v^e' 


2ia 


sina 


4 Case n = 2 

4.1 The sub-case of 8 real roots 


We will use the reduction formula (1.6) applied to the three integrals, assuming as before 1 < a < 6 , we get: 

r.2 


(a, 6) := 


X" 


(4.1) 

(4.2) 

(4.3) 


and reply the path done in theorems |3.2[ |3.5| and |3.7| obtaining further formulae which provide tt in theorems 

[TH|0| e[131 


^0 ^(x 2 -a 2 )(x 2 -A)(:r2_52)(^2_^) 

sR 2 ^°'’°'Ha,b) := [ . ^ dx 

Ji f 

/ OO 

- V(^ 


/(^2_a2)(^2_ l)(^2_52)(^2_ 1) 


a;2-a2)(x2- A)(a:2-62)(a:2- T) 


2_L 


11 


















































Interval [0,1/b] 

Computing (a, &) we get: 

Theorem 4.1. If 1 < a < & then 


TT = 


263 

p(3) / 

3.111 

2 ’ 2 ’ 2 ’ 2 


1 A 

" l 

2 

"O 

a262 ’ 64 j 


K 


'b{o? - 1)' 


1 1 \ V ~ ^ I ® ~ 1) / +1 


K 


a (&2 + 1) 


Proof. Applying the reduction (1.61 one finds, (we omit at all details for shortness: see theorem 3.2), 

'b{a^-l)\ 


si?' 


[ 0 . 1 / 6 ] 


1 


{a,b) = - 


1 


2 U2 _ 1 


K 


b (g^ + 1) 
' (52 — 1 ) / 6^ + 1 "" I a + 1) 


K 


Thesis (4.4) follows by lemma 2.1 which provides: 


8-^2 


[ 0 . 1 / 6 ] 




3.111 
2 ’ 2 ’ 2 ’ 2 


1 1 

b‘^ ’ 


(4.4) 


(4.5) 


(4.6) 

□ 


As for corollary 3.3 for a —> 1 by (4.4) we get a formula for tt with only one parameter: 


Corollary 4.2. If & > 1, then: 


25^(62 - 1 ) 


K 


Vl+6:^ ) 


1 + 6^ / 1, ^ 

64 + (1-62)fJ 2^2 


(4.7) 


jp 

’ b^ 


Interval [1 /a, a] 

Computing 8 i? 2 ^^°’°^ (a, &) we get: 
Theorem 4.3. If 1 < a < 6 , then 


a (h^ — 1) 


K 




7'(3) 

D 


1 • _ 1 1 1 
2 ’ 2 ’ 2 ’ 2 


1 ^4 (a‘^-l)&^ a^-1 


(4.8) 


Proof. For computing 8 i? 2 ^^"^’'*^ (a, 5) we meet the same feature seen for the proof of 3.5 where the employ of 
Cauchy-Schlomilch transformation requires to split the integration domain.We are then led to: 


8??^^/“’“' (a, b) = ^ 


6^ — 1 y a (6^ — 1) 
On the other side, by formula (2.4 d) proved in lemma [2)l]ve have: 


84^/“’“'(a,fe) = 


•nab 


2\/b'^ — aP'y/aP'b'^ — 1 


A3) 

D 


1 • _ 1 1 1 
2 ’ 2 ’ 2 ’ 2 


1-a^ 


- 1 ) - 1 


a'^ — b'^ ’ 0^6^ — 1 


Thesis (4.8) follows by comparing (4.9) and (4.10). 


(4.9) 


(4.10) 


□ 


Remark 4.4. Let us note that, even with reference to our previous papers iniiniEn] how in ( |4.8[ ) is the first 
appearance of a Lauricella function whose parameters bi are not all equal. Then the variables order cannot be 
changed. 


Small changes to theorem 3.7 proof, allow, starting by integral 8 i? 2 '’°°^(aj b), to obtain the third formula to 


n when the polynomial (1.8) has 8 real and different roots. 
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Interval [6, oo) 

Theorem 4.5. If 1 < a < 6 then: 


TT = 


,_ (6' + l)K 

V(62 - a2) {a?h^ - 1) 




l)K 


\a{h^ + l)) 


ab'^y/b^ — 1 

f ( 3 ) / 

1.111 

2 ’ 2 ’ 2 ’ 2 

1 1 \ 


" 1 

1 

1 — 64 ’ 1 — afb'^ ’ — b^ J 


(4.11) 


Proof. By applying (1.3) to the integral 
creasing, we get: 


{a,b), taking into account the Cauchy-Schlomilch is now in- 




1 


2 \ 62 _ 1 


b ^ - 1) 




a (62 — 1) y 62 + 1 \ a (62 + 1) 


(4.12) 


On the other side, through formula ( 2.5|3 ) (proved in lemma 2.1) also we get: 


Rp°^\a,b) = 


irab^ 

2^{b^~-iy{b^ - a2) (a262 -1) 


/ 1.111 


1 


1 


1 — 6^ ’ 1 — a262 ' — b'^ 


As usual, thesis (4.11) follows by comparing (4.12) and (4.13). 


(4.13) 


□ 


In this case again we can take the limit for a —>■ I"*" getting by (4.111 a formula to tt with only one parameter 


like the formula (3.20) shown in 3.8 


Corollary 4.6. If 6 > 1 we have 


2(1 + +k(^) 


1 -64 + 26V&^ - 1Fi f ^’1’ ^ 


(4.14) 


1 - 62 ’ 1 - 64 


Proof. If in (4.11) we take the limit for a —>■ I"*", since one of two elliptic integrals becomes of module zero and 
then degenerates in a elementary one, see equation (3.10), we obtain: 




1 - 64 + 26- 1 f 2>2^2>2 


1 


1 


1 


1 - 62 ’ 1 - 62 ’ 1 - 64 


Thesis (4.14) follows by identity (3.9) already recalled for in the proof of corollary 3.20 


(4.15) 


□ 


4.2 Sub-case with 4 real and 2 pairs of conjugate complex roots 


If the P(x)’s structure is (3.21), we know that P{x) > 0 for a: G [0,1/6] U [6, oo) and then we will consider the 
integrals: 


(a, 6) = 


0 — 5^) — 62) (^2 — e“2*a) 


/ OO 


dx 


: dx 


— 5^) (x^ — b^) (a;2 — e“2*a) _ e2“) 

Through the usual approach of a two-fold evaluation (lemma |2.2[) we will find two further formula to tt. 


(4.16) 

(4.17) 
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Interval [0,1/b] 


We will start from (4.16). 

Theorem 4.7. If 6 > 1 then 

2b‘^ (^{1 + b'^) K 


_ 2b sin a _ 

Vb^ — 2b^ cos 2a+l 


) - V54-262cos2o+lK(^|f^)) 


(l + 6^) \/b^ — 262 cos 2a + 1F 
Proof. By lemma [2^ we have 


' 3.111 

(3) / 2’ 2’ 2’ 2 

^ ' 2 


(4.18) 


64 ’ 62 ’ 62 




64 ’ 62 ’ 62 


(4.16i) 


But the transformation (1.3) drives as usually, 6) to two integrals both linked to the entry 3.152-12 

p. 277 of [g: 




26 sin a ^ 6 ® 

2 V \/64 — 262 pQg 2a -I-1 V \/64 — 262 cos 2a + 1 / 1-1-62 _j_ ^2 


K 


( | 4 T 6 ] 3 ) 

□ 


Thesis (4.18) follows by (4.103) and (4.103). 

Interval [b, 00 ) 

Like (4.17), evaluating 4 i? 2 ^’°°^(a,6), equation (4.17) we get 
Theorem 4.8. If a, 6 S K, 6 > 1 then 

y^b'^ — ly/1 — 262 pQg 2a + 64 


/ 


/ 1.111 

^2pg) 2^Y>2 


IT 


26 sin a 


1 1 1 \ 

1 — 64 ’ 1 — 62e“2»a ’ 1 _ 62e2®“ / 3 


•\/l — 262 cos 2a + 64 
Vl — 262 cos 2a + 64 




26 cos a 
1 - 1-62 


1-^62 


Proof. Thesis follows because lemma [2?^ establishes 

b^n 


4Rl^'°°\a,b) = 


/ i- 1 1 1 

_p(3) / 2’ 2’ 2’ 2 

2-\/&4 _ x-y/l — 2 cos 2a 62 -|- 64 ^ \ 1 


1 — 64 ’ 1 — 62e“2*a ’ 1 — 62f 


But through the reduction scheme induced by transformation (1.3), it is also true that: 

6 / 26 cos a' 


= ) (j- 




K 


-1-62 \ 1 -I- 62 / Vl — 2 cos 2a 62 -|- 64 V Vl — 2 cos 2a 62 -|- 64 


26 sin a 


then (4.19) follows by comparison of (4.173) with (4.17 3c). 


(4.19) 

( |4T7l 3) 

jilTl oc) 

□ 


4.3 Sub-case: 4 pairs of coniugate comple roots 

Interval [6, 00 ) 

In such a case we will get a result similar to|3.12[ The path is now done more quick being easy to check that: 


; da; = 


da: 


Jo a/( 1 — 2a;2 cos 2a -I- a;4) (1 — 2a;2 cos 26 -|- a;4) Jo a/( 1 — 2a;2 cos 2a + x'^) (1 — 2x‘^ cos 26 -|- a;4) 

(4.20) 

namely, in our notation: 

(a, &) = {a,b)- dF^ ) 

The integral o7?g°’°°^(a, 6) at the right side of (4.20) has been computed during the proof of theorem 
formula (3.30i). Going into details, we have: 


3.12 
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Theorem 4.9. Let a, 5 S K be such that sin a > sin 6. Then: 


K 


\/ sin^ a — sin^ b 


sin a 


2 sin a / 2 ’ 2 ’ 2 ’ 2 > 2 


l + e2“,l + e-^“,l + e^*",l + e 


^ — 2ia 


2ib 


^—2ib 


Proof. By lemma [273| taking s = 2 we obtain 

/ 1. 1 1 1 1 
4 


p[U,OOJ/ 7\ _ TT'V 

ori2 (o, o) — ^ r p 


1 + 6^“, 1 + e-2“, 1 + 1 + 


On the other side by (4.20) and (3.30d it follows that 


o4°’“^(a,6) = 


1 f \/ sin^ a — sin^ b 


2 sin a 


K 


sma 


(4.21) 


(4.22) 


(4.23) 


Then thesis (4.21) follows by comparing!4.22) with (4.23). 


□ 


Remark 4.10. Both thesis (3.31) and (4.21) of theorems 3.12 and 4.9 differ only in parameter “a” of Lauricella’s 
function. In other words theorems 3.12 and |4.9| drive to the conclusion that, due to the peculiar choice of the 
arguments (l + e^“, 1 + 1 + 1 + e~^) both Lauricella functions will coincide. 


5 Case n = 4 


In the case with exponent 4 we use the relationship (1.7) and this will require a greater circumspection due to 


the fact that the Cauchy-Schlomich transformation applied to Roberts’s integrals, leads to a difference of two 
divergent integrals. In addition, the exponent 4 at numerator does not allow any integration on unbounded 
intervals. 


5.1 Sub-case: 8 real roots 

Here we will face only the integrals: 

:= f ^ dx (5.1) 

•^0 ^(:,2_„2)(^2_^)(^2_^,2)(^2_^) 

8 i ?4 ^“’“'(O) ^) := / , ^ da; (5.2) 


being ^I^l'°°\a,b) = +oo. 


Interval [0,1/b] 

Theorem 5.1. If 1 < a < 6 then: 


8i?f'/''(a,6) 


1 


jb^-b^ + l) ni-a2)b \ 

( 1 - 62 ) yail-b-^)) 


(6^ + 62 + l) / (a2 + i) \ 

6(1 + 62) l^a(62 + l)y 


-(1 


6^) E 


\a{l-b^)) 


(1 + 62) ^ 


/(1^\ 

l^a(l + 62) j 


(5.3) 
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Proof. First we will carry out the integration on [e, 1/5]; after the transformation (1.31, we will take the limit 
for e —>■ O'*'. Then we get: 


5) = lim 

e->0+ , 


= lim 

£->■0 + 


^ (x2 - a?){x‘^ - i)(a;2 - }P){x‘^ - i) 

l-u^ 

A+i 


dx 


-drt 


fa + r) — 


1 f-. 

2A-i 


l + v"^ 


(a - i) i) 

We must take the limits ought to both integrals of kind u> a> fi 


\ 

dt! 


(5.4) 


G{u,a,l3) = f 

J a 


a/ {a^ — x‘^){l3'^ — x'^) 


dx 


appearing in (5.41 are divergent for x —>■ +oo and in (5.4) one shall evaluate their difference. Now, by the entry 
3.153-9 of [13] we have: 

G{u, a, (3) = a 


u (3 


(3 a 


F arcsin 7 ;, — — E \ arcsin —, — 


u [3 


/?’ a 




Recalling the integration formula (3.6), thesis (5.3) follows. 


□ 


Using again lemma [2T[ through the usual double evaluation, we find a formula providing tt where complete 
second kind elliptic integrals appear too. 

Theorem 5.2. If 1 < a < b then: 


8 4 IK -|- E 

TT = -b -;;;- 

3 F 


ove 


and 


finally: 


-b"^+ 1) ^ ( {a^ - 1) 5^ (54 + b‘^ + l) ^f{a^ + l)b 


52 - 1 


:(52 - 1) 


1-^52 


K 


a (1 -I- 52 ) 


(5.5) 

(jK^) 


E = (1 - 5^) £; 


2^ ^ I , A , r2^ ^ / (a2 + l)5 


a (52 — 1 ) 


-f (1 + 5^) E 


a (1 -I- 52 ) 


/ 5 . 1 1 1 

F = FgM 


1 1 


52 ’ 0252 ’ 5"^ 


(1^) 

@) 


Proof. Through the lemma 2.1 for computing (a, 5) we get: 

o / 5 . 1 1 1 

n[0.1/&]/ (3) / 2’ 2’ 2’ 2 

8^4 (0,5)- 144 I 


1 1 


52 ’ 0252 ’ 5^ 


Thesis (5.5) follows by comparison (5.3) to (5.6) and getting tt. 


(5.6) 


□ 


In such a case too, it makes sense to take the limit for a —>■ I"*" and obtain a formula holding only one 
parameter where the Lauricella function degenerates in an Appell one. 

Theorem 5.3. If 5 > 1, then 


85' 


d {b^-l)K[j^) + {l-b^){l + b^yE[j^) 


l + b^ 


456-h3(l-52)Fi 


5.1 1 

2 ’ ’ 2 


(5.7) 


62 ) 64 
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Interval [1 /a, a] 


By the transformation (1.3) we compute the integral (5.2) and obtain: 


Theorem 5.4. Ifl < a < b then: 




1 


6(62 _ 1) 


(1 - 6^ + b'^)K 


Proof. By (1.3) we get: 




i(62 - 1) 


l + v'^ 


- {V - lyE 


- 1 ) 


z(62 - 1) 


■ dt! 




Thesis (5.8) follows from (5.8 d), with the help of entries 3.153-5 p. 277 and 3.152-7 p. 266 of 


By lemma (2.1) and theorem 5.4 we obtain: 


Theorem 5.5. If 1 < a < 6 then: 

(1 - (gfeg) - {b^ - l^E 


62 (62 - 1) 


1. _ 3 1 1 

,(3) / 2> 2’ 2’ 2 

^ 1 


1 _ g^-l 

^ ’ a2-62 , a^b^-i 


5.2 Sub-case: 4 real and 2 pairs of conjugate complex roots 

Interval [0,1/b] 


In such a case P{x) has the shape (3.21), then, via the (1.3), we are going to study the integral: 

rr „,4 


•^0 y (a;2 — (a;2 — 62) (1 — 2cos2aa:2-I-a;^) 

By means of our usual technicalities, we get the reduction: 

Theorem 5.6. Se6>leaGKsi ha: 


dx 


1-^62 


(^4 _ ^2 ^ 

-I — - = IK 2 (a, 6) — \/l — 262 cos 2a -I- 6^ £ 2 ( 0 ,6) 

vl — 262 cos2a -I- 6"^ . 


where 


- / i.\ „/26cosa\ , ,, ^f2bcosa 

^i(o,6) = , Ei(a, 6) = £; 


l-t62 


l2ia,b) = K 


26 sin a 


v/l — 262 cos 2a + b^ 


1 -h62 
£ 2 ( 0 , b) = E 


26 sin a 


Vl — 262 cos 2a -I- b'^ 


By thesis (5.11) and by lemma 2.2 form. (2.6), we conclude: 


Theorem 5.7. If 6 > 1 and a G M we get: 

86^ 


3Fg) ( 

5.1 1 1 

2 ’ 2 ’ 2 ’ 2 

1 e-2>“ 1 


f,4 , (,2 j &2 1 

\ 

3 


(1 + 62) El (a, 6) - 


1 + 62 


^ +1) ^ _ y^i_ 262 cos 2a+ 64 £ 2 ( 0 ,6) 


•\/l — 262 (,Qg 2a + 6"^ 


(5.8) 


pb) 


□ 


(5.9) 


(5.10) 


(5.11) 


(5.Ill) 


(IpTja) 


(5.12) 
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6 Conclusions 


In this article based upon a Roberts’s reduction approach of hyperelliptic integrals to elliptic ones, and on 
the simultaneous multivariable hypergeometric series evaluation of them, several identities have been obtained 
expressing tt in terms of special values of elliptic, hypergeometric and Gamma functions. We acted selecting the 
values n = 0, 2, 4 in equation (1.21 and to the roots of reciprocal equation P{x) = 0, where P{x) is provided by 


(1.81. Accordingly, we obtained formul® where tt can be constructed trough only one parameter as in formulae 
(3.11), (3.20), (4.7), or through two parameters like (4.4). Let us note the 7r-formula ( 3.31| ) holds two parameters 
and the imaginary unit too. 

In any case, such results will be a useful tool in order to check the routines which one can build for 
the practical computations of Lauricella’s functions we too met frequently in our researches on Mechanics 
or Elasticity, as in PBIHllg]. 

These relationship are, as far as we are concerned, all unpublished and undoubtably unknown not only to 
any human being but also to any computer algebra systems like Mathematica(g). 

For example, by implementing our (3.11) in Mathematica(g we get: 


ln[1]:= pp[b_] : = 


2*2 (-1 + *2) EllipticK[--i^] 


(l + *2) +2 (-l + *2) AppellFl[i, 1, 1, ^]) 


in[2]:= FullSimplify [pp [b] ] 

2 b2 ( - 1 + b2 ) EllipticK [ ] 

Out[2]= -^- - -- 

-bMl + b2) +2 (-l + b^) AppellFl[i, 1, b, 1, 


Figure 2: Formula (3.11) 


Mathematical is not able to reduce the above expression to tt. Nevertheless proceeding numerically we see 
that the results is confirmed within the range of numerical turbulence. 

in[3]:= Max [Table [tt - pp [b] , {b, 1.001, 20, .05}]] 

Out[3]= 4.12061 X 10-“ 

in[4]:= Plot [Abs [tt - pp [b] ] , {b, 1.001, 20}, PlotStyle-» Black] 


2.x 10*''^ 



Figure 3: Numerical check of formula (3.11) 
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